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Summary 

This paper presents the results of applying the methodology of singular perturbations and time 
scales (SPATS) to the control of digital flight systems. A block-diagonalization method is described that 
decouples a full-order, two-time-scale (slow and fast) discrete control system into reduced-order, slow 
and fast subsystems. Basic properties and numerical aspects of the method are discussed. This study 
reveals an interesting fact that singularly perturbed discrete systems can be viewed as two-time-scale 
systems. The closed-loop optimal control of the full-oider, two- time- scale system involves the solution 
of a full-order, algebraic, matrix Riccati equation. Alternatively, by using the block-diagonalization 
method described in this paper, the full-order system s decomposed into reduced-order, slow and fast 
subsystems. The closed-loop optimal control of the subsystems requires the solution of only reduced- 
order, algebraic, matrix Riccati equations. A composite, closed-loop, suboptimal control system is 
constructed as the sum of the slow and fast, optimal fe idback controls. The application of this method 
to an aircraft model shows close agreement between tbe exact solution and the decoupled (composite) 
solution. The main advantage of the method is the considerable reduction in the overall computational 
requirements for the evaluation of optimal guidance and control laws. The significance of the result is 
that it can be used for real-time onboard simulation. This paper also contains a brief survey of digital 
flight systems. 

1. Introduction 

The dynamics of many systems is described by high-order differential equations. Frequently, the 
presence of small parameters such as time constants, masses, and moments of inertia is the source for 
the increased order of the system. A system in which the suppression of a small parameter is responsible 
for the degeneration of the dimension of the system U called a singularly perturbed system. Such a 
system possesses widely separated clusters of eigenvalues exhibiting “slow” and “fast” phenomena. 
The high dimensionality coupled with the two-time-scale (slow and fast) behavior makes the system 
computationally “stiff,” with the result that extensive numerical routines are required. 

The theory of singular perturbations and time scales (SPATS) in continuous control systems has 
reached a level of maturity (refs. 1 to 3) and has been successfully applied to aerospace problems (refs. 4 
to 15). On the other hand, the subject of SPATS in dig tal flight control systems has not received much 
attention so far. 

The purpose of this paper is to present a methodo ogy of SPATS to discrete control systems with 
an application to an aircraft model. The basic ides, underlying this method is to apply a block- 
diagonalization procedure to decouple a full-order, two Time-scale system into low-order, slow and fast 
subsystems. This procedure offers a considerable redu tion in the overall computational requirements 
for developing optimal guidance and control laws. The significance of the result is that it has possible 
applications for real-time onboard implementation. Also, this paper contains a brief survey of digital 
flight systems. The results presented in this paper are m expanded version of those recently published 
by the authors in references 16 to 19. 

The organization of the paper is as follows. Section 2 describes the concept of digital flight control 
and the previous research in this field. The ideas of SPATS in continuous and discrete control systems 
are given briefly in section 3. Section 4 presents a bloc.<-diagonalization procedure to decouple a high- 
order, two-time-scale system into low-order, slow and fast subsystems. An aircraft example is provided 
to illustrate the method. In section 5, the focus is on tbe time-scale synthesis of digital, optimal control 
systems. The closed-loop, linear, optimal regulation of the original two- time- scale system requires the 
solution of a high-order, algebraic, matrix Riccati equation. Alternatively, a composite, closed-loop, 
suboptimal control is constructed as the sum of slow and fast, optimal feedback controls, which require 
the solution of only low-order, algebraic, matrix Riccati equations. Numerical results obtained for an 
aircraft example show excellent agreement between th;‘ original (optimal) solution and the composite 
(suboptimal) solution which is computationally simple \ 
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Subscripts: 

c 

/ 

9 

s 
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input matrix 

input matrix for subsystem (i = 1,2) 

matrix of order (m X n) 

matrix of order (n x m) 

eigenvalue of ( ) 

feedback matrix 

identity matrix 

performance index 

discrete interval 

order of state, z 

order of state, x 

positive, definite symmetric matrix of Riccati equation 
eigenvalue 

positive, semidefinite symmetric matrix of order (n-j-m) x (n+m) 
positive, definite symmetric matrix of order (r x r) 
order of control, u 
time 

control vector of order r 

state vector of order n 

output vector 

state vector of order m 

small, positive perturbation parameter 

composite mode 
fast mode 
correction mode 
slow mode 
zero order 


Superscripts: 

T transpose 

(0) zero-order solution 

Abbreviations: 

ADC analog-to-digital converter 

DAC digital-to-analog converter 

max maximum 

min minimum 

A dot over a symbol indicates a derivative with respect to time. 
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2. Background 

In this section, the idea of digital flight control systems is discussed briefly and their advantages 
are enumerated. The previous research in digital flight ontrol is summarized. 

A digital control system uses a digital computer to mplement the control logic. The development 
of reliable, faster, and inexpensive microcomputers has recently aroused considerable interest in digital 
control systems (refs. 20 to 24). The first commercial digital computers in the 1950’s were not fast 
enough or small enough to be placed on space vehicles. The NASA Apollo program spurred the 
development of smaller and faster computers for the digital control of boosters and spacecraft in 
the 1960’s. This technology was transferred to aircraft by NASA in the 1970’s. The invention of 
microprocessors in the 1970’s made digital computers very small, feist, and inexpensive. Many military 
and civilian aircraft now have digital control systems (ref. 25). Digital control, which has made possible 
more accurate and sophisticated autopilot logic, promises to be the focal point in most improvements in 
navigation, guidance, and control in the future. A digital flight control system contains analog-to-digital 
and digital- to-analog converters (ADC and DAC, respe< tively) as shown in figure 1. 



Figure 1. Concept of digital flight control. 

With adequate redundancy incorporated into the design, digital flight control systems ensure 
adequate flight safety. The present capabilities for incorporating integrated circuits into lightweight, 
low-cost minicomputers and microcomputers make digital implementation of modern flight control 
systems especially attractive (refs. 25 to 28). 

Another feature of digital implementation is the potential for the synthesis of complex control 
systems that involve high-order nonlinearities and that utilize time-sharing for multiple-loop control. 
One such complex control structure is an adaptive system that is capable of online adjustment of the 
control parameters in response to changing flight charac eristics (ref. 29). 

The advantages of digital flight control systems over t heir continuous counterparts are summarized 
(ref. 30) as follows: 

(1) Use of complex control strategies 

(2) Easy implementability 

(3) Flexibility in adding or changing functions 

(4) Repeatability of performance 

(5) System integration and hardware economy 

(6) Reduction in cost, size, weight, and power dissipation 

In the past, there have been many applications of digital control theory for flight systems. The first 
application of digital technology to flight control was digital implementation of basic analog autopilot 
functions (refs. 31 and 32). Modal control theory has been applied to the design of digital flight control 
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systems such as pitch-attitude control systems and roll/yaw control systems for a short takeoff and 
landing (STOL) aircraft (refs. 33 to 36). An algorithm was presented for obtaining a reduced-order 
or simplified nonlinear model of the F-8 aircraft (ref. 37). A method for the synthesis of a nonlinear, 
automatic, flight control system involving an F-8 aircraft has been developed (refs. 35 and 38). A 
residue-measure criterion has been employed as an efficient method for model reduction in the design 
and analysis of the digital flight control system of NASA’s Space Shuttle (ref. 39). Discrete design 
methods have been compared at various sample rates for several different flight modes for the terminal 
configured vehicle (TCV) at the Langley Research Center, which is a modified Boeing 737 aircraft 
(ref. 40). Research has been conducted by NASA into some advanced control laws for the F-8 digital 
fly- by-wire (DFBW) program (refs. 29 and 41). 

The analysis and control of large systems has been a formidable task, not only because of the high 
order of the systems but also because the majority of these systems possess interacting phenomena of 
different (slow and fast) speeds. The simultaneous presence of slow and fast dynamics makes the system 
computationally “stiff,” with the result that extensive numerical routines are required. A stiff system, 
having a two-time-scale character, need not necessarily be in the singularly perturbed structure. 


3. Singular Perturbations and Time Scales 

In this section, the definitions of singular perturbations and time scales in continuous and discrete 
control systems are given briefly. Broadly speaking, a continuous (discrete) system described by a 
differential (difference) equation containing a small parameter is called a singularly perturbed system 
if the order of the system is reduced by neglecting the small parameter. Such a system possesses widely 
separated clusters of eigenvalues and exhibits a two-time-scale (slow and fast) character. 


Continuous Control Systems 

Consider briefly the idea of SPATS in continuous control systems and its application to aerospace 
problems. Consider a linear time- invariant system such as 


x — Anx + A 12 Z + Biu 

C-K 

II 

3 

II 

0 

(la) 

ei = A 21 X + A 22 Z + B 211 

(* (t = 0 ) = z(o)) 

(lb) 


where x and z are n- and m-dimensional state vectors, respectively, and u is an r-dimensional control 
vector. The matrices A„ and B; are of appropriate dimensions, and £ is a small, positive perturbation 
parameter. The (n -j- m) high-order system (eqs. (1)) is in the singularly perturbed form in the sense 
that by making e = 0 in equations (1), the resulting system is of reduced order n only, and consequently 
the degenerate system does not satisfy the initial condition z (0) . That is, 

i = A u x (0) + Ai 2 z (0) + Biu x (0) (« = 0) = x(0) (2a) 

0 = A 2 ix(°) + A 22 z (0) + B 2 u (»<°> (t = 0) £ z(0)j (2b) 

Here, let us note that equation (2a) is a reduced-order system of order n only, with equation (2b) being 
an algebraic equation, and the initial condition z(0) is lost in the process of degeneration. The singular 
perturbation theory retrieves the lost initial condition. The block diagrams of the full, high-order, 
or perturbed system (eqs. (1)) and the degenerate, low-order, or unperturbed system (eqs. (2)) are 
shown in figure 2. Here we assume that the input u is independent of e. Otherwise, u becomes u ■ . 

Also, x and z are predominantly slow and fast state vectors, respectively. In a general two-time-scale 
(slow and fast) system, the singular perturbation parameter e does not appear explicitly as shown in 
equations (1). The essence of the theory of SPATS is to take advantage of the order reduction associated 
with the degeneration and to decouple the original system (eqs. (1)) into slow and fast subsystems of 
reduced order for further analysis and synthesis. 
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(a) Original system. 



(b) Degenerate ; ystem. 


Figure 2. Singularly perturbe i continuous system. 

The theory of SPATS has been successfully applied to aerospace problems described by continuous 
control systems (refs. 4 to 15). The digital flight control systems with singular perturbation (or two- 
time-scale) character have not received much attention so far. 

Discrete Control Systems 

For the purpose of synthesizing digital controllers, he continuous model (described by differential 
equations) of a dynamic system is transferred into a discrete model (described by difference equations) 
where the controls are held constant between the sampling intervals (refs. 34, 42, and 43). 

The general form for a linear, shift-invariant, singularly perturbed discrete system is given by 
(refs. 44 to 47) 

x(/c+ 1) = An x (*0 + £ 1- 1 A 12 z(fc) + Bi u(fc) (3a) 

e 2i z(k+ 1) = £ J A 2 ix(fc) + e A 22 z(fc) +£ j B 2 u (k) (3b) 

with 0 < i < 1 and 0 < j < 1 where x(/c) and z k) are slow and fast state vectors of n- and 
m-dimensions, respectively, u(fc) is an r-dimensional control vector, e is the singular perturbation 
parameter, and A and B t are matrices of appropriate dimensionality. The initial-value problem with 
x(fc = 0) = x(0) and z(A: = 0) — z(0) is formulated. 
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The three limiting cases of equations (3) result in the following models: 


(1) For the C-model (i = 0; j = 0), 

x(A:+l) — An x(fc) 4 eA \2 z (k) 4 Bi u (k) (4a) 

z(fc+ 1) = A 21 x(fc) 4 eA 22 z (k) 4 B 2 u (k) (4b) 

where the small parameter e appears in the column of the system matrix. 

(2) For the R-model ( i = 0; j = 1), 

x(fc-bl) = An x(fc) -f A 12 z (k) 4 Bi u(fc) (5a) 

z(fc4l) = £A 2 i x(fc) 4 eA .22 z(fc) 4 eB 2 u (k) (5b) 

where the small parameter e appears in the row of the system matrix. 

(3) For the D-model (i = 1; j = 1), 

x(fc4l) = An x(fc) 4 A 12 z(A:) 4 Bi u(Ar) (6a) 

ez(k + l) = A 21 x(Ar) 4 A 22 z(fc) 4 B 2 u(k) (6b) 


where the small parameter e is positioned in a fashion identical to that of the continuous systems 
described by differential equations. 

For the present discussion, it is enough to consider the system (eqs. (4)) as an initial- value problem 
with x(/c = 0) = x(0) and z (k = 0) = z(0). The suppression of the small parameter e in equations (4) 
results in the degenerate system 

x( 0 )(Hl) = A n x(°>(i) + Biu(fc) (x<°>(Jk = 0) =x(0)) (7a) 

* ( °)(fc+l) = A 2 ix(°)(fc) + B 2 u(fc) (zW(fc = 0)/z(0)) (7b) 

Here equation (7a) is a difference equation in x(°)(/c) of order n only, whereas equation (7b) is an 
algebraic equation. It means that once x^(k) is solved from equation (7a), z(°)(fc) is automatically 
fixed by equation (7b) and zi°i(0) is not, in general, equal to z(0). Thus, the suppression of the small 
parameter e in equations (4) leads to a low-order system (eqs. (7)) with a consequent loss of the initial 
condition z(0). Hence, by definition, equations (4) are in the singularly perturbed form. The original 
system (eqs. (4)) and the degenerate system (eqs. (7)) are shown in figure 3. Here it is also assumed 
that u is independent of e\ otherwise u (k) becomes u W(fc) when e is made equal to zero. The theory 
of SPATS in discrete systems is of recent origin (refs. 44 to 55). Attempts have been made to apply 
the technique to digital flight control systems, but these attempts have been limited mainly to a class 
of digital control of continuous systems (refs. 56 to 61). 

In this section, the concepts of SPATS in continuous and discrete control systems have been 
discussed. In the next section, a method for decoupling a two- time- sc ale system is presented. 

4. Time-Scale Analysis of Digital Control Systems 

In this section, attention is focused on time-scale analysis (TSA) of digital control systems. This 
is followed by a description of a block-diagonalization procedure to develop a method to decouple a 
discrete system into slow and fast subsystems. The basic properties and numerical aspects of these 
systems are discussed. An interesting fact is provided showing that singularly perturbed discrete 
systems can be viewed as two-time-scale systems. Finally, an aircraft example is given to illustrate the 
method (ref. 16). 
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(a) Original system. 



(b) Degenerate system. 


Figure 3. Singularly perturb i d discrete system. 


Two-Time-Scale System 

Consider a linear, shift-invariant, two-time-scale discrete system 


x(fc + l) 


An A12 

'x(fc)' 

4- 

Bi 

_ z(fc-l-l) _ 


A21 A22. 

z(fc) 

I 

. B 2 . 


u(fc) 


(8) 


with initial conditions x(fc = 0) = x(0) and z(fc = 0) - z(0). It is assumed that the system (eq. (8)) 
is asymptotically stable and that its eigenspectrum consists of a cluster of n large eigenvalues and a 
cluster of m small eigenvalues. Let the eigenvalues of equation (8) be arranged as 


l| > |Psl| > • ■ • > |Psn| > P/ll > • • ■ > |P/ml 


(9) 


If the condition 

e = |p/il/|psn « 1 (!°) 

is satisfied, then the system (eq. (8)) possesses a two- time-scale (slow and fast) property (refs. 44 
and 50). In other words, the stable discrete system (eq. (8)) is said to exhibit a two-time-scale behavior 
if the largest absolute eigenvalue of the fast eigenspectrum is much smaller than (i.e., has wide sepa- 
ration with) the smallest absolute eigenvalue of the slow eigenspectrum. In the literature, one usually 
finds stiffness being stamped on a system in which the ratio of eigenvalues is given by equation (10). This 
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ratio appears explicitly as the small parameter e in the case of a singularly perturbed system. It is 
also noted that the slow modes are generated by large eigenvalues and the fast modes are produced by 
small eigenvalues. 

In order to block-diagonalize the system (eq. ( 8 )) or separate it into slow and fast subsystems, a 
two-stage linear transformation is used. In the first step, the block A 2 i is removed to make equation (8) 
an upper-triangular matrix by using the transformation 


where the (m x n) matrix D is a real 


Zf(k)=z(k) + D x(fc) 
root of the Riccati-type algebraic equation 


A 22 D — DA|i + DAi 2 D — A 2 i = 0 

Substituting equation ( 11 ) into ( 8 ) results in 


where 


x(k+ 1) 


A s Aj2 

' x(*) 

1 

B, 

Zf(k+1) _ 


0 Ay 

Z/(fc) 

-+■ 

[B/J 


u(k) 


A s = An — A 12 D Ay = A 22 + DA 12 

By — DBj + B 2 

In the second stage, the transformation 


(ii) 


( 12 ) 


(13) 


x a (fc) = x(fc) + Ezf(k) ( 14 ) 

is applied to equation (13) and the (n X m) matrix E is chosen as the solution of the Riccati-type 
algebraic equation 


EAy — A S E + A\2 — 0 (15a) 

E(A 22 + DAi 2 ) - (An — Ai 2 D)E + A 12 = 0 (15b) 

Then, equation (13) reduces to 


where 


x s (A;-)-l) 


A fl 0 ] 

'x s {k)' 

i 

B s 

_Zf (k+l) _ 


< 

o 


T 

[ByJ 


B s = (I 5 + ED)Bi + EB 2 


(16) 


Now the system (eq. (16)) is in the desired decoupled form, as shown in figure 4 . Note that the slow 
and fast variables x s (/c) and zy(fc) are related with the original state variables x(/c) and z (A;) by means 
of the transformations given in equations (11) and (14), which are combined to form 


Xg(fc) 


' (Is + ED) E 

'x(fc)' 

z f \k)_ 


D r /J 

z(k) 


and from which 


'x(fc)' 


Ce 

1 

H 

Xs(k)' 

.»(*) . 


[-D (Iy + DE)J 

Z/(fc) 


(17) 


(18) 


Note that the special feature of these transformations is that equation (18) needs no inversions of the 
matrices, and hence the simplified computation is derived. 
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Figure 4. Decoupled slow ind fast subsystems. 

The condition given in equation (10) can be written as 

e = max|e(Ay)|/mia|e(A 5 )| « 1 ( 19 ) 

where 

e ( A /) = {Pf 1> P/m} 

e(A s ) = {p s 1, • ■ , Psn} 

The iterative solution of equation (12) is given as (ref. 44) 

D J+1 = (A 22 D t + D, A 12 Dj - A 21 ) A”/ (20) 

with an initial value of Do = — A 2 iA 11 1 . Similarly fo" equations (15), we have 

Ei+i = Afj 1 (E t A 22 + EjD Ai 2 + Ai 2 DE t - + Ai 2 ) (21) 

with an initial value of E 0 = AJ’/Aij. Substituting these initial values of D 0 and E 0 into equation (16) 
gives 


A a o = A-ll + A 12 A 21 A i / 

A /0 = A 22 - A 2 iA 11 1 Ai 2 

(22a) 

B s o = Bi - A^j 1 Ai 2 A 22 A 11 1 Bj 

+ A^ 1 A 12 B 2 

(22b) 

Byo = B 2 — A 2 jAjj Bj 


(22c) 

X s o(fc+l) _ A s0 0 

_Zf 0 (k + l)\ 1 0 A /0. 


(23) 


Thus, it is seen that a two-time-scale system is decor pled into slow and fast subsystems. 
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Singularly Perturbed Systems as Two-Time-Scale Systems 

We consider the singularly perturbed system (eqs. (4)) and repeat it here for convenience. Thus, 


x(fc+ 1) 


An 

eA i2 

'x(Ar)' 

1 

Bi 

z(fc + l) _ 


. A 21 

eA 2 2_ 

z(k) 

+ 

B 2 


(24) 


with initial conditions x(/c = 0) = x(0) and z (k - 0) = z(0). 

Now, substituting the transformation of equation (17) into (24) and replacing E by eE allows the 
decoupled subsystem to become 


+ 

no 

X 

1 


A s 

0 

x s ( k)' 

1 

B s 

L z /(*+i). 


0 

£ A f 


+ 

[B/J 


where 

As = Aj i — £Ai 2 D A/ = A 22 + DA 12 

B s = (I s + eED)Bi + eEB 2 B / = DB 1 +B 2 
and D and E are the solutions of the Riccati-type algebraic equations 


(25) 


eA 22D — DAn + sDAi 2 D — A 21 = 0 
eE(A22 + DA 12 ) - (An ~ eAi2D)E + A 12 = 0 

whose iterative solutions start with the initial values of D 0 = — A 2 iAj _ 1 1 and E 0 = A _1 Ai 2 . For a 
sufficiently small e, the zero-order approximate solutions are n - 


(26) 

(27) 


Aso = An 
B s o = Bi 

so that the decoupled subsystem becomes 


A/o = A22 - A 2 iA 11 1 Ai2 

B/o = B2 — A 2 iA 1 " 1 1 Bi 


x s0 (fc + l) 


'An 

0 

X s o(*)' 


Bi 

Zyo(fc+l) 


0 

L 

£ A/o. 

Zfo(k) 

+ 

-B/o. 


Xso {k = 0) = x(0) 

z f0 (k = 0) = z(0) - A 2i Aj“ 1 1 x(0) 


(28a) 

(28b) 


(29) 


Now the system (eq. (29)) is obtained in an alternative way from equation (24) if the singular 
perturbation method (SPM) is used. The degenerate (slow) subsystem is obtained by making e = 0 in 
equation (24). That is, 


x(°)(fc+l) = Aiix( 0 )(fc) + B,u(°)(fc) (x(°>Ofc = 0) = x(0)) (30a) 

z^°)(/c+l) = A 2 i + B 2 u^°\k) (z( 0 )(A: = 0) z(0)) (30b) 

Here, it is noted that z (k) has lost its initial condition z(0) in the process of degeneration. In order 
to recover this lost initial condition, a correction subsystem is used. The transformations between the 
original and correction variables are 


Xg(k) = x(k)/e M z g(k) = z{k)/e k 

u g (k) = u {k)/e^ 1 

Substituting equations (31) into (24) allows the transformed system to become 

£x 9 (fc + l) = An x 9 (Ar) -(- Ai 2 z g (k) -f Bi u g (k) 
z g (k + l) — A 2 i x 9 (fc) + A 22 z g (k) + B 2 u 9 (A:) 


(31a) 

(31b) 
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(32a) 

(32b) 



The zero-order approximation (e = 0) of equations (32) becomes 


0 = A u x' 0) (/c) + A - 2 z5 0) (fc)+Bi u^\k) (33a) 

4 0) (fc + l) = A 21 x^ 0) (fc) + A 2i zi 0) (fc) + B 2 4 0) (A:) (33b) 

Rewriting equations (33) produces 

Xg°^ (k) = -Aj'j 1 A12 Z;°^(A:) + Bi u 9 °^(fc)j (34a) 

z^ 0) (fc + l) = A 90 Z 9 °V) + B 9 o 4 0) (A:) (34b) 


where 


A 9 o = A 22 - A2 i A 11 I Ai2 
B 9 o = B2 - A21 Ajj'Bi 

The total solution consists of the sum of the outer solution given by equations (30) and the correction 
solution of equations (34). For the present, to simplify the presentation, u(k) and its associated 
functions are omitted. Then (ref. 47), 

x(/c) = x^(fc) + e** 1 x^(fc) (35a) 

z(k) = zW(k) + e k z { ° ) (k) (35b) 

The' current interest is only zero-order approximations, i.e., the slow part x(°)(fc) ofx(fc) and the fast 
part e k z { g\k) of z (k). Thus, from equations (30) and (34) (omitting input functions), 



x(°)(/c+l) = An X* ] \k) 

(36a) 


Zg ^(fc+1) = AgO Z'g ^ {k) 

(36b) 

or 

z) 0) (fc+l) = eA 9 o4 0) (fc) 

(36c) 

where 

Z r°\k) = £ k Zg^lk) 

(37a) 


x (°)(fc = 0) = x(0) 

z (°)(/c^o) = zJ 0) (o; 

(37b) 


= z(0)-z(°)(0) 

= z(0) - A 2 i Aj'j 1 x(0) 

(37c) 


By comparing equations (29) and (36), it is seen that to a zero-order approximation, the degenerate 
and correction subsystems of the singular perturbat on method are the same as the slow and fast 
subsystems of the time-scale analysis. This is a result similar to that in continuous systems where the 
singular perturbation and time-scale approaches give identical results to a zero-order approximation. 
The method is now demonstrated by an aircraft example. 

Application to an Aircraft 

The methodology described in the previous sections is illustrated by the application to an aircraft. 
The longitudinal dynamic equations of an aircraft model are obtained based on a rigid-body assumption 
(ref 62) Also, the angle of attack is assumed to be small. The linearized model possesses a two-time- 
scale property in the sense that pitch angle, velocity and altitude are the slow variables and angle 
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of attack and pitch rate are the fast variables, corresponding to phugoid and short-period modes, 
respectively. For digital implementation, a zero-order hold is used to obtain the discrete model. The 
aircraft under consideration is a twin-engine general aviation aircraft. The discrete model is given by 
(ref. 63) 


y(fc + l) = Ay(fc) + Bu(fc) 


(38) 


where 


y(fc) = 



A - 

A n 

A 12 

D _ 

sr 

z(fc)_ 

tx. — 

A 21 

A 2 2_ 

D — 

B 2 . 


and 


Xl(/c) 

velocity, ft/sec 

x 2 {k) 

pitch angle, deg 

x 3 (fc) 

altitude, ft 

zi{k) 

angle of attack, deg 

*2 (k) 

pitch angular velocity, deg/; 

ui(k) 

elevator deflection, deg 

u 2 (k) 

flap deflection, deg 

u 3 (k) 

throttle position, deg 


B 



' 0.923701 

-0.308096 

0 

0.053043 

-0.090367- 


0.039705 

0.995525 

0 

-0.107454 

0.588883 


0.087127 

1.899490 

1.0 

-0.635270 

0.394015 


-0.035537 

0.010123 

0 

0.007748 

0.137407 


. 0.069562 

-0.012706 

0 

-0.097108 

0.287411. 


r 0.042825 

-0.000395 

-0.154048- 




-0.484628 

-0.515359 

-0.002237 



= 

-0.161525 

-0.067522 

-0.005257 




-0.202010 

-0.289303 

0.005061 




.-0.806770 

-0.852161 

-0.006353. 




The states xj (k) and x 3 (/c) are scaled down by a factor of 100 to facilitate easy representation. The 
eigenvalues of the discrete (original) system are given by 1.0, 0.962103 ± 0.175342j, 0.217297, and 
0.072882, suggesting that there are three slow modes (xi(fc), X 2 (k), and x 3 (Ar)) and two fast modes 
(xi(fc) and z 2 (k)) with the eigenvalue separation ratio given as |p/i|/|p s 3 | = 0.222196. The results are 
now summarized below for two cases in which the two-time-scale approach is used. 


Case I — Zero iterations. Here the initial values Dq and Eo are themselves used in obtaining 
the decoupled subsystems without resorting to any iterative solutions of the Riccati-type algebraic 
equations (26) and (27). Thus, 


D n = 


Eo = 


A s o = 


0.038398 

0.001715 

0 

-0.074861 

-0.010405 

0 

0.021141 

0.098166 ' 


-0.108780 

0.587615 


-0.430486 

-0.730705 


0.914899 

-0.309128 

0 

0.087915 

1.001837 

0 

0.141016 

1.904675 1.0 
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A /o 


0.009600 0 134947 
-0.099961 0 288048 


B s o = 


-0.040446 

-0.935733 

0.513256 


- 0.089654 -0.153556' 

- 0.981356 0.000913 

0.676144 -0.008690 


B/o 


-0.201197 

-0.804933 


-0.290202 

-0.846769 


-0.000858 

0.005203 


The eigenvalues of the slow subsystem are 1.0 and 0.9 >8368 ± 0.159020/, and the eigenvalues of the 
fast subsystem are 0.225595 and 0.072053 with an eigei value separation ratio of 0.232222. 


Case II Five iterations. Here the values of D5 and E5, obtained after five iterations of the Riccati 
equations (26) and (27), are used for the decoupled subsystems. The following matrices are obtained: 


D S 


0.023371 -0.006771 0 

-0.110005 -0.028>61 0 


E 5 


-0.026374 0.215647 

-0.178397 0.741358 

-0.144412 -1.487336 


A s 5 — 


0.912520 

0.106996 

0.145317 


—0.310327 

1.011676 

1.906477 


0 

1.0 

1.0 


A /5 


0.009715 0.131307 

-0.099864 0.280473 


B S 5 = 


-0.123959 

-1.041054 

1.053315 


-0.173428 

-1.093552 

1.219171 


-0.151789 

0.005406 

-0.021322 


B 


/ 5 


-0.197728 —0.285322 0.001476 

-0.797591 -0.837347 0.010657 


The eigenvalues of the slow subsystem are 1.0 and 0.962098 ± 0.175345/, and the eigenvalues of the 
fast subsystem are 0.217307 and 0.072882 with an eigenvalue separation ratio of 0.222207. 

The summary of eigenvalues for the original system an d the decoupled subsystems is shown in table I. 
The eigenvalue separation ratio (consisting of the maximum of absolute values of fast eigenvalues and 
the minimum of absolute values of slow eigenvalues) is s mwn in table II. 

The responses of the original system and the decoupled subsystems are shown in figures 5 to 9 
for case I with zero iterations and in figures 10 to 14 for case II with five iterations. From both the 
figures and the tables, it is clear that the decoupled subsystems are close to the original system with 
considerable reduction in total computational requirements. 

In this section it has been shown that a full-order, two-time-scale system can be decoupled into 
reduced-order, slow and fast subsystems. The technique has been demonstrated by applying it to an 
aircraft model. The next section considers the decoupling of optimal discrete control systems. 
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Table I. Summary of Eigenvalues Table II. Summary of Eigenvalue Separation Ratio 



Bap ‘0|6ue ipRd 



14 


Figure 5. Comparison of solutions for velocity for case I. Figure 6. Comparison of solutions for pitch angle for case I 
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5. Time-Scale Synthesis of Digital Optimal ( ontrol Systems 

In this section a two-time-scale discrete control system is considered. The closed-loop, optimal, 
linear, quadratic regulator for the system requires the solution of a full-order, algebraic Riccati equation. 
Alternatively, the original system is decomposed into reduced- order, slow and fast subsystems. The 
closed-loop optimal control of the subsystems requires the solution of two algebraic Riccati equations of 
lower order than that required for the full-order system. A composite, closed-loop suboptimal control is 
formed from the sum of the slow and fast feedback opt mal controls. The main advantage of the method 
is a considerable reduction in the overall computational requirements for the closed-loop optimal control 
of digital flight systems (refs. 18 and 64 to 67). 

For the sake of convenience, let us repeat here tiie stable, linear, shift-invariant discrete control 
system given as equation (8): 


x(/c + l) 


A 11 a 12 

’x(fc)' 

— 1— 

b r 

z(fc + l) 


. A 21 A 22 . 

z(fc) 

I 

b 2 . 


u (k) 


where the various state variables and the matrices have been previously described. The performance 
index to be minimized is 

J = if^[y T (k)Qy(n + u T (k)Ru(k)} (39) 

1 k=0 


where 


y T (k)= [x T (I), z T (k) 


and 


Qll Ql2 

Q 12 Q22 . 


is a positive, semidefinite symmetric matrix of orde (n + m) x (n + m). The matrix R denotes a 
positive, definite symmetric matrix of order (r X r). 

The closed-loop optimal control is given by (ref. 68) 


u(Jfc) = -R _1 B r P (l + BR _1 B T P) 1 Ay(fc) 


(40) 


where P, of order (n + m) x (n+m), is the positive, definite symmetric solution of the algebraic, matrix 
Riccati equation 

P = A r P (l + BR“ B r p) _1 A + Q (41) 

where 

a _ r a h A t2 

a 21 A 22 

The closed-loop optimal system is given by 

y(fc+l) = (A- BF)y(fc) (42) 


B = 


Bi 

B 2 


where 

F = R _1 B r P (l + BR~ 1 B r P)~ A 

Instead of directly tackling the original regulator problem described by equations (8) and (39), 
one should appropriately decompose it into two regulator problems (approaches) for slow and fast 
subsystems. For this, it is necessary to separate the original performance index into the sum of two 
performance indices for slow and fast subsystems. 


17 



Separation of Original Performance Index 


Approach 1 . In approach 1 the original performance index (eq. (39)) has to be represented as the sum 
of the performance indices of the slow and fast subsystems. Using the transformation in equation (18) 
between the original state variables (x(fc) and z(fc)) and the subsystem variables (x s (Ar) and z f(k)) in 
equation (39) and using u s (k) — u f(k) = u (k) as shown in figure 4 results in 


where 


.j ^ 

J = - ^ [xj (k) Q ss x s (fc) + Zf {k) Qff z f{k) + uj (fc) R a u 5 (fc) 
Z 0 

+ xf(fc)Q a/ zy(fc) + zJ{k)Q fs x s (A:) +uj(k)R f uy(A:)j 


Qss — [Is, D 


Qii Q12 


I S 

. Q12 Q22 . 


-D 


(43) 


Q// = [E T , -(I / + DE) 5 


Q11 Q12 


E 

. Ql2 Q22 . 


-(I/ + DE) 


Q sf = [Is, D' 


-Q11 Q12 

E 

Ql2 — Q22 . 

(1/ + DE) 


-Qll Q12 


V 

Ql2 -Q22. 


D 


Q/s = [E r , (I/ + DE) 5 


Note that Qf s = Qjj and R s = = R/2. 

Since J has to be represented as the sum of J s and Jy, it is necessary to neglect Q s f and Q f s 
Then, 


i 

J s = - ^ [xf (k) Q ss x s (fc) + uj (fc) R a u s (fc)| (44) 

^ 0 
1 OC 

jj = 2 Ys [ z /( fc )Q// z /( fc ) + u /( fc ) R / “/(*)] ( 45 ) 

As Q s f and Q f s have been neglected, this certainly introduces an error in that J will not be equal to 
the sum of J s and Jj. To take this into account, it is necessary to readjust Q ss and Qff (ref. 65). 

Approach 2. Since the present method is simply a design or synthesis problem, first select the 
performance indices of the subsystems and then formulate the total performance index of the original 
system. Thus, define 


1 oo 

J 's = [ x ^ ( fc ) Qss *•(*) + (*) R s u-(*)l (46) 

^ 0 
1 OO 

J f = 2^2 [*/(*) Q//*/( fc ) + u /( fc ) R / u /( fc )] ( 47 ) 
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Using the transformation in equation (17) between the subsystem variables 
original system variables x(k) and z (k) produces 


x s (/c) and Zf{k) and the 


1 ^ 

= 2 X] [y r ( fc ) Q' y( fc ) + uT (k) r' u(fc)j (48) 


where 


Q' = 



q; 2 

Q22. 


(r' = 2R' S = 2R' / j 


Q'll = (Is + ED) r Q^(I s +ED) +D r Q' // D 
Q'i 2 = (Is + ED) r Q' s E + B T Q' ff 
Q22 = d t qUe + q} / 


Thus, in equations (46) to (48), first select Q' s , Q' ff , R' , and R}; and then by using D and E, one 

obtains Q' and R'. Here it is possible to decouple J J s , and J f exactly without any approximation. 

The original J is dependent on D and E, the decoupling matrices, which may not be of practical 
advantage. 


Optimal Control of Subsystems 

Using the transformation in equation (17) allows i he original system in equation ( 8 ) to decompose 
into slow and fast subsystems as 


x s (/c + 1 ) = A s x s (i:) + B s u s (fc) ( 49 ) 

z/(k+l) = AfZf(k) + B f u f (k) (50) 

Now these slow and fast subsystems are optimized with respect to their corresponding performance 
indices given in equations (44) and (45), respectively. The slow regulator problem consists of the slow 
subsystem of equation (49) and the performance ind.-x of equation ( 44 ). The fast regulator problem 

consists of the fast subsystem of equation (50) and the performance index of equation ( 45 ) For 

convenience, let Q ss = Q s and Q// = Q f . 

The optimal feedback control of the slow subsystem of equation (49) is given by 

u s (fc) = -R-JbJp, (i s +B,R s - 1 Bj’P s ) _ 1 A s x s (A:) (51) 

where P s , of order (n x n), is a positive, definite symmetric solution of the reduced-order, algebraic 
Riccati equation 

P s = A S P S (l s + B S R S J BjP s j A s + Q s (52) 

Similarly, the optimal feedback control of the fast subsystem of equation (50) becomes 

u f (k) = -R^BjP; (iy + B j f )~ l A f z f (k) (53) 

where P y, of order (m x m), is a positive, definite symmetric solution of the reduced-order, algebraic 
Riccati equation 

P / = A / P / (l/ + ByRy ’BjPy) Ay + Qy ( 54 ) 
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( 55 ) 


The control laws of equations (51) and (53) are rewritten, respectively, as 

Ug(fc) = -F s x s (fc) 


u f {k) = -FfZ f {k) 


(56) 


where 


F s = R~ l BjPs (l s + B s R s Ag 

F f = RJ l B T f P f (l f + B f R-J 1 B T f Pf) l A f 


It is noted that the control laws in equations (51) and (53) are optimal only with respec t to t he slow 
and fast subsystems of equations (49) and (50). It is, however, computationally simpler to deterame 
these control laws than the optimal control law of equation (40) of the original system (refs. 51 and 52). 

Composite Control 

The composite control is formulated as the sum of the slow and fast feedback controls given by 
equations (51) and (53). That is, 


u c (fc) = u s (fc) + U/(fc) 

= - [F,x s (fc) + Ff zy(fc)] 


(57) 


Applying the transformation of equation (17) between the slow and fast variables and the original 
variables to equation (57) produces 

(58) 


uc(fc) = - [F C s x(fc) + F cf z (fc)] = — F c y(fc) 


where 


(59) 


Fcs = F s (Ig + ED)+F / D 

F c/ = F S E + F ; 

F c = [F C s, F c/ ] 

Using the composite control of equation (58) in the original system in equation (8) results in 

y c (fc+ 1) = (A - BF C ) y c {k) 

as shown in figure 15. It is known that minimizing the original performance index of equation (39) with 
respect to the 8 composite system of equation (59) results in the suboptimal performance index (refs. 65 

and 69), 

Jc = \ y r (o) P c y(o) ( 6 °) 

where P c is the positive, definite symmetric solution of the discrete Liapunov equation 

P c = (A - BF c ) T P c (A - BF C ) + Q + FjRF c (61) 
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(a) Optimal feedl ack control. 

x c (K) 
z c (k) 


(b) Composite feedback control. 

Figure 15. Original and composite systems. 

Aircraft Example 

Consider the same aircraft model as previously given. The eigenvalues of the discrete (original) 
model are 1.0, 0.962103 ± 0.175342j, 0.217297, and 0.072882, an indication that there are three slow 
modes (x\ (fc), X 2 (fc), and X 3 (k)) and two fast modes (z\{k) and 22 (h)) with an eigenvalue separation 
ratio of 0.222196. The performance measures in eq lation (39) are taken '‘as Qn = Q 22 = I and 
R* = R / = 1. 

Using the method described in the previous section gives the results summarized below. The 
positive, definite matrix P of equation (41) is 

■ 4.822583 0.427518 0.709917 -0.272287 0.059216- 

0.427518 11.654952 3.732418 -2.665579 3.042291 

P = 0.709917 3.732418 2.690923 -1.096067 0.882337 

-0.272287 -2.665579 -1.096067 1.742674 -0.689301 

. 0.059216 3.042291 0.882337 -0.689301 2.015052. 

Using the above P, the optimal control of equation (40) is obtained with F as 

‘ -0.021105 -1.136159 -0.208959 0.260873 -0.507318 - 

F — -0.110915 -0.700095 -0.101750 0.142104 -0.410844 

-0.629305 -0.025190 -0.100651 0.033007 0.003334 
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The closed-loop optimal system of equation (42) has the eigenvalues 0.843545, 0.258790 ± 0.326649;, 
and 0.06457 ± 0.047911j. 

For the slow and fast subsystems, the corresponding values are 


Ps 


4.621454 

1.051044 

0.945863 


1.051044 

13.073895 

4.463674 


0.945863 

4.463674 

2.710497 


and 


p / = 


1.204794 

0.24435 


0.244345 

1.914426 


For the composite control, the feedback matrix F c , obtained from equation (58), is 


F c 


-0.010457 

-1.138702 

-0.270112 

0.262312 

-0.507373 

-0.110238 

-0.699239 

-0.101678 

0.142112 

-0.410528 

-0.654236 

-0.018681 

-0.100062 

0.031065 

0.004602 


and the positive, definite matrix P c , obtained from the discrete Liapunov equation (61), is 


Pc 


■ 5.045733 

0.338736 

0.671223 

-0.205383 

-0.013507- 

0.338736 

11.734399 

3.769142 

-2.716177 

3.103791 

0.671223 

3.769142 

2.714962 

-1.128617 

0.918364 

-0.205383 

-2.716177 

-1.128617 

1.790493 

-0.744196 

.-0.013507 

3.103791 

0.918364 

-0.744196 

2.083318. 


The corresponding closed-loop, suboptimal composite system has eigenvalues 0.840132, 0.259351 ± 
0.334107;, and 0.063358 ± 0.04656;. The performance indices of the original optimal system and the 
composite suboptimal system are obtained as 612.1122 and 616.1734, respectively, with a percentage 
error of 0.66347. The central processing time on the CDC CYBER computer required for the original 
system is 5.42 sec, whereas that required for the composite system is only 0.44 sec. Throughout the 
paper, the computations are performed by using the linear-quadratic-regulator (LQR) design package 
ORACLS and commercial mathematical library (ref. 70). 

The responses for the various states and controls for the exact (optimal) system and the composite 
(suboptimal) system are shown in figures 16 to 23. These responses are obtained for the case when the 
aircraft trajectories are to be regulated to the equilibrium conditions of x\ e = 190.66 ft/sec, X 2 e — 0, 
X 3 e = 2000 ft, and z\ e = Z2e = 0 with equilibrium controls. 

The above results clearly show an excellent agreement between the exact system and the composite 
system. Once again it is noted that the composite control is obtained from the lower-order, slow and 
fast subsystems with a considerable reduction in the overall computation (refs. 16 to 18). 


Concluding Remarks 

This paper describes a time-scale method for digital control systems with an application to an 
aircraft. The idea of digital flight control systems was briefly discussed and their advantages were 
enumerated. The previous research in digital flight control was summarized. The concepts of singular 
perturbations and time scales (SPATS) in continuous and discrete control systems were introduced. 

Attention was focused on the time-scale analysis of digital control systems. A block-diagonalization 
procedure was described for developing a method to decouple a discrete system into slow and fast 
subsystems. Basic properties and numerical aspects of the method were discussed. It was shown that 
singularly perturbed discrete systems can be viewed as two-time-scale systems. Finally, an aircraft 
example was given to illustrate the method. 

The important problem of optimal control of two-time-scale discrete systems was discussed. The 
closed-loop, optimal, linear, quadratic regulator for the system requires the solution of a full-order, 
algebraic, matrix Riccati equation. Alternatively, the original system was decomposed into reduced- 
order slow and fast subsystems. The closed-loop optimal control of the subsystems requires the solution 
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Figure 18. Comparison of solutions for altitude. Figure 19. Comparison of solutions for angle of attack. 
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Time, sec Time, sec 

Figure 22. Comparison of solutions for flap deflection. Figure 23. Comparison of solutions for throttle position. 


of two algebraic Riccati equations of an order lover than that required for the full-order system. 
A composite, closed-loop, suboptimal control has been formed from the sum of the slow and fast 
optimal feedback controls. Again, numerical results obtained for an aircraft model showed very close 
agreement between the exact (optimal) solution and the composite (suboptimal) solution which is 
computationally simpler. The main advantage of the method is a considerable reduction in the overall 
computational requirements for obtaining the closed- loop, optimal control laws of digital flight systems. 
The significance of the methodology is that it can be used for generating optimal guidance and control 
strategies for real-time onboard simulation. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
September 15, 1988 
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